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1. Introduction 

The European Twisted Mass Collaboration (ETMC) is currently performing large scale sim- 
ulations with A'y = 2 + 1 + 1 flavors of dynamical quarks using Wilson twisted mass lattice QCD 

Particular problems are caused by the non-degenerate strange/charm quark doublet, since be- 
sides the usual twisted mass parity breaking the strange and charm quark numbers are not con- 
served. The latter amounts to contamination of correlation functions by intermediate states with 
wrong flavor quantum numbers. In particular, correlation functions of charmed mesons and baryons 
obtain small contributions from similar strange systems, which are significantly lighter and, there- 
fore, inevitably dominate at large temporal separations. In this paper we discuss these problems 
in the context of the kaon and D meson masses, which are important quantities already at the 
stage of generating gauge field configurations, when tuning the strange and charm quark masses to 
their physical values. We present three methods to overcome these problems and demonstrate their 
consistency. This work is a summary of a more detailed recent paper 

2. Simulation setup 

The gauge action is the Iwasaki action The fermion action is Wilson twisted mass, 

SF.nght[X^'\x^'\U] = a'''£x^'\x)(^DMmo) + 1^75^3) X^'\x) (2.1) 

SF,he..y[X^'\x^'\U] = a4£^W(x)(Dw(mo) + /McT75Ti+T3At5)z^''W (2.2) 

for the light degenerate up/down doublet and the heavy non-degenerate strange/charm doublet 
[^] respectively, where Dw denotes the standard Wilson Dirac operator. K = 1/ (2mo + 8) is tuned 
to maximal twist by requiring m^^,f"^ = 0, which guarantees automatic ^(a) improvement for 
physical quantities, e.g. the here considered kaon and D meson masses. We also refer to [^, where 
this Nf = 2+1 + 1 twisted mass setup has been pioneered. 

All results presented in the following correspond to computations on 1042 gauge field con- 
figurations from ensemble B35.32 characterized by gauge coupling j8 = 1.95, lattice exten- 
sion L? xT = 32^ X 64 and bare untwisted and twisted quark masses K = 0.161240, pL = 0.0035, 
/ifj = 0.135 and jig = 0.170. The corresponding lattice spacing is a 0.078 fm, the pion mass 
mps 318MeV. 

3. Quantum numbers, physical and twisted basis meson creation operators 

In Wilson twisted mass lattice QCD parity is not a symmetry and the heavy flavors cannot be 
diagonalized - both symmetries are broken at (^(a). Consequently, instead of the four QCD heavy- 
light meson sectors labeled by heavy flavor and parity, {s,—), {s, +), (c, — ) and (c, +), there is only 
a single combined heavy-light meson sector {s/c, — /+) in twisted mass lattice QCD. In contrast to 
QCD, where the D meson is the lightest state in the (c, — ) sector, it is a highly excited state in the 
combined {s/c, — /+) sector of twisted mass lattice QCD. This in turn causes severe problems for 
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the computation of mo, since the determination of excited states is inherently difficult, when using 
lattice methods. 

To create kaons and D mesons a suitable set of operators is 



(3.1) 



In the continuum the twisted basis quark fields X and their physical basis counterparts y are related 
by the twist rotation 



■(«) 
(d) 



y(") 

exp(/75T3«//2) I 



(.0 



X 



(■v) 



exp(/75TiW/,/2y 1 ^(c) 



• (3.2) 



0)/ and (Oh denote the light and heavy twist angles, which are 7r/2 at maximal twist. At finite lattice 
spacing the procedure is more complicated, since renormalization factors have to be included. The 



heavy-light meson creation operators ( |3.1[ ) and the corresponding renormalized operators in the 
physical basis transform into each other via 



/ + 



ClCh -SlSk -SlCh -ClSh 



\ 



-SiSh +ClCh -ClSh -SlCh 
+SlCh +ClSh +ClCk -SlSh 
\ +ClSh +SlCh -SlSh +ClCh J 



+zsx^''h^''^ 
\ -zsx^'h^'^ J 



(3.3) 



where Cx = co?,{(Ox/2), = sm{(Ox/2) and Zp and Z5 are operator dependent renormalization con- 
stants. 

The starting point for the three analysis methods presented in the following sections ^ to ^ are 
the 4x4 correlation matrices 



Cjk{t) = {a\ffj{t){ffm) 1^) 



t, 



(3.4) 



of spatially extended, i.e. APE and Gaussian smeared versions of twisted basis heavy-light meson 
creation operators (p|). The smearing parameters have been optimized by minimizing effective 
masses at small temporal separations (cf. [HI] for details). 



4. Method 1: solving a generalized eigenvalue problem 

One possibility to determine niK and mj) is to solve the generalized eigenvalue problem 



Cju{t)^^f(t,t^ = Cjk{t^)X^"Ktfyy;'{t,to), 



(4.1) 



cf. e.g. [P] and references therein. From the eigenvalues A'") one then computes four effective 
masses <ffective by solving 



AW(?,?o) 



eM-m'-^ecti^e i^Mt) + exp(-'W^ffective ^o) - 0) 



(n) 



Xi"){t + l,^o) exp(-m(l,,,(^^o)(^ + 1)) +exp(-m(l,,,(f,fo)(r - {t + 1))) 



, (4.2) 
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Figure 1: left: effective masses obtained from the 4x4 correlation matrix (3.4) by solving the generalized 
eigenvalue problem (4. 1 ) and (4.2) and corresponding plateau fits; right: zoomed effective mass of the lowest 
state, which has been identified as the kaon. 



where T is the temporal extension of the lattice and the exponentials e{—m^^^^^:^^^{t,to){T — t)) and 
^(~"^effective(^^o)(7' ~ {t + 1))) take care of effects due to the temporal periodicity of the lattice. 
Heavy-Ught meson masses are finally determined by fitting constants to effective mass plateaus at 
temporal separations t ^ \. Results are shown in Figure |l|. 

It is well known that for sufficiently large temporal separations t the generalized eigenvalue 



problem ( [4.1| ) and ( |4.2[ ) yields the lowest states in the sector considered, i.e. in our case the four 
lowest states in the combined {s/c,—/+) sector. The D meson (m{D) 1868 MeV), however, is 
not among them. Lighter states include 

• the kaon and its radial excitations, 

m{K) w 496MeV, m(i«:(1460)) = 1400 MeV - 1460 MeV, ... 

• parity partners of the kaon 

m{K*{SOO)) = 672(40)MeV, m{K*{U30)) = 1425(50)MeV, ... 

• multi particle states 

m{K + K), m{K + l x k), ... 



At first glance it seems that the 4x4 correlation matrix (3.4) is not sufficient to determine md, but 
that one needs a significantly larger correlation matrix, which is able to resolve all states below the 
D meson. Note, however, that in the continuum an exact diagonalization of Cju is possible yielding 
one correlator for each of the four sectors {s,—), {s,+), (c, — ), (c,+). Hence the generalized 
eigenvalue problem would not yield the four lowest masses of the {s/c, — / +) sector, but tuk, '«(i,+)> 
niD and At finite lattice spacings such a diagonalization is, of course, only approximately 

possible. However, discretization artefacts, which are responsible for that, only appear at 0'{a) 
and are thus expected to be small. Therefore, at not too large temporal separations one of the four 
effective masses should be dominated by the D meson and, consequently, provide an estimate for 
niD. 

One can check that this is indeed the case by twist rotating the eigenvectors v'"' to the pseudo 



physical basis, which is defined by (3.3) with Zp = Zs = 1. Then one can read off the approximate 
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flavor and parity content of each of the four states corresponding to the four effective mass plateaus. 
As shown in Figure the lowest state is dominated by the physical basis operator y^^^ Ys 
therefore, interpreted as the kaon, while the second excited state is dominated by i/A^^^^ 75 y/'' \ i.e. 
corresponds to the D meson. 
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Figure 2: left: squared pseudo physical basis eigenvector components of the lowest state 'p, which is 
interpreted as the kaon; right: squared pseudo physical basis eigenvector components of the second excited 
state I v^-^' p, which is interpreted as the D meson. 



5. Method 2: fitting exponentials 

An alternative approach to determine mk and mj) is to perform a minimizing fit of 

£ (4"V4"' (exp ( - m„?) + exp ( - m„(r -?))), (5.1) 

i.e. of A'^ exponentials to the computed correlation matrix Cjk in a suitably chosen window of tem- 
poral separations t. Notice that T is the temporal extension of the lattice and the second exponential 
exp(— m„(r — t)) takes care of effects due to the lattice temporal periodicity. The masses m„ can 

in) 

be interpreted by analyzing the prefactors Uj , which is very similar to what has been explained 
in more detail in the previous section for the eigenvector components provided by the generalized 
eigenvalue problem. 

As before we find that the lowest state is a kaon, while the second excited state is dominated 
by the physical basis operator xj/^'^^ 75 i/A^^' and, hence, should correspond to the D meson. 



6. Method 3: heavy flavor/parity restoration 



Our third approach is based on the twist rotation of heavy-light meson creation operators (3.3). 
In a first step we express the correlation matrix Cjk in the physical basis in terms of the twist angles 
ft)/ and (Oh and the ratio of renormalization factors Zp and Z5: 

CP^y^^c-^'\t;(o,,(OH,Zp/Zs) = 

= ^{o)i,0}h)diag{Zp,Zp,Zs,Zs)C{t)diag{Zp,Zp,Zs,Zs)^''{(Oi,0),,), (6.1) 
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where the matrix ^ has been defined in (p3|). We then determine CO/, (Oh and Zp/Zs by requiring 

cf''''''\t;couCOH,Zp/Zs) = 0. (6.2) 

At finite lattice spacing and small temporal separations t this cannot be achieved exactly, because 
of iff (a) heavy flavor and parity breaking effects and the presence of excited states. However, at 
sufficiently large t, where only the kaon survives, the condition (62) can be realized. It amounts to 
removing any kaon contribution from the diagonal correlators (^P^y**"^^'^^ j ^ (^^^ _y 

Finally we analyze the diagonal correlators (jP^y^^^'^^'^ individually. There is one correlator for 
each of the four sectors {s, — ), {s, +), (c, — ) and (c, +). The effective mass plateaus corresponding 
to the (s,—) and the (c,— ) diagonal correlator yield the heavy-light meson masses and nio, 
respectively. 



7. Conclusions and outlook 

Results for mk and mo obtained with our three methods agree within statistical and system- 
atic errors, see Table |l] and Figure |3[ For a detailed discussion, of how statistical and systematic 
errors have been determined, we refer to [Q]. We are able to determine mk in a rigorous way 
with rather high statistical precision (statistical error <0.4%). On the other hand all three methods 
require assumptions, when computing nio, which amount to a systematical error being involved. 
Nevertheless, the combined statistical and systematical error for nij) is <2.5%. 





method 1 


method 2 


method 3 


rriK 
mo 


0.2184(3) 
0.829(8) 


0.2177(8) 
0.835(20) 


0.2184(3) 
0.823(15) 



Table 1: comparison of kaon and D meson masses determined with the three methods presented. 



mass of the K 



0.2195 
0.219 



method 1 (generalized eigenvalue problem) < 
method 2 (fitting exponentials) ' 
method 3 (heavy flavor/parity restoration) ■ 



2 

method 



mass of the D 



0.87 
0.86 
0.85 
0.84 
0.83 
0.82 
0.81 



method 1 (generalized eigenvalue problem) ■ 
method 2 (fitting exponentials) 
method 3 (heavy flavor/parity restoration) - 



2 

method 



Figure 3: comparison of kaon and D meson masses determined with the three methods presented. 



Being able to determine mk and md is very important for tuning the strange and charm quark 
masses to their physical values. For precision charm physics, however, we intend to use a mixed 
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action Osterwalder-Seiler setup [|T0[]. First steps in this direction are currently under way and have 
been reported during this conference 
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